Collisional decoherence of a tracer particle moving in one dimension 
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We study decoherence of the external degree of freedom of a tracer particle moving in a one 
dimensional dilute Boltzmann gas. We find that phase averaging is the dominant decoherence effect, 
rather than information exchange between tracer and gas particles. While a coherent superposition 
of two wave packets with different mean positions quickly turns into a mixed state, it is demonstrated 
that such superpositions of different momenta are robust to phase averaging, until the two wave 
packets acquire a different position due to the different velocity of each wave packet. 
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I. INTRODUCTION 

The transition from quantum mechanics to classical 
physics is one of the most debated problems in the his- 
tory of modern physics. In particular, the question arises 
why one can not observe macroscopic objects in a super- 
position of spatial distinct locations, despite the fact that 
all objects are made up of microscopic particles which 
indeed can be observed in such position superposition 
states. Several conceptually very different solution to this 
problem were proposed, as e.g. the theory of spontaneous 
localization pQ, which modifies the Schrodinger equation 
by adding an incoherent part. A less drastic approach 
is the theory of environmentally induced decoherence [2] . 
This assumes that the combined state of system and en- 
vironment evolves according to Schrodinger's equation, 
but if only the system density operator is observed, it 
seems as if the coupling to the environment destroys the 
quantum feature that a system can be in a superposition 
of several distinct states, a process known as decoherence. 

During the last two decades, there has been increas- 
ing interest in the engineering of large quantum systems, 
e.g. for quantum information processing. A major lim- 
itation to these efforts is posed by their fragileness to 
decoherence. Therefore, a detailed understanding of dif- 
ferent decoherence processes is no more just an academic 
problem, but necessary for future quantum technologies. 

A paradigm of environmentally induced decoherence 
is collisional decoherence, where the system of interest 
is a tracer particle, possibly macroscopic in size, which 
experiences random collisions with particles of a ther- 
mal reservoir. The colliding particles can be molecules, 
much as in Brownian motion, but one could also think of 
photons of the cosmic background radiation. Several au- 
thors put forward increasingly complicated master equa- 
tions for a tracer particle in a thermal gas, first for an 
infinitely heavy tracer particle [3-6 , and later for a tracer 
particle with finite mass [7HTU]. The latter were used to 
study collisional decoherence by applying quantum tra- 
jectory methods [TTJ [12]. However, the validity of the 
single collision calculations used in the derivations of the 
respective master equations for tracer particle with finite 
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FIG. 1: (color online) (a): Undisturbed double slit experiment 
with screen in the far field (not to scale), (b): A colliding 
particle measures through which slit the tracer particle passes, 
removing the coherence pattern at the screen, (c): A colliding 
particle transfers momentum to the tracer particle, shifting 
the phase of the coherence pattern. If this shift is random, the 
coherences can not be observed anymore. One of the results 
of this paper is that process (c) is responsible for decoherence, 
rather than process (b). 



mass was recently questioned [13] [14] , and a consensus is 
still missing. 

The mechanisms by which an environment can destroy 
measurable superpositions can roughly be divided into 
two categories. First, the environmental state can get en- 
tangled with the system state. This effectively delocalizes 
the relative phase of any superposition state of the sys- 
tem into the combined state of system and environment. 
After tracing out the environmental degrees of freedom, 
this leads to a reductions of the coherences of the sys- 
tem. We say, the environment measures the system (see 
figure [I] (b)). In the second process, sometimes called 
phase averaging, the interaction with the environment 
changes the relative phase of the superposition state. If 
this phase change is random and different for each run 
of an experiment, then the coherence of the superposi- 
tion state can no longer be observed (see figure [l](c)). In 
some sense, phase averaging does not fundamentally de- 
stroy the coherences, but rather makes the relative phase 



unpredictable and therefore unmeasurable. 

One should mention that neither mechanisms com- 
pletely solves the problem of the quantum-classical tran- 
sition because the measurement problem remains [2]. 
Nevertheless, both mechanisms can successfully describe 
the observed lack of coherences of macroscopic objects. 
However, it is the first decoherence mechanism which is 
mostly cited in connection with the quantum-classical 
transition, possibly because is appears to destroy coher- 
ences more fundamentally compared to phase averaging. 
This view is especially maintained within the topic of 
collisional decoherence [3j HJ [15] . 

Although a colliding gas particle certainly carries away 
some information about the state of the tracer particle, 
we show in this article that the decoherence due to this 
information exchange is negligible compared to phase av- 
eraging. The latter arises because a collision adds a rel- 
ative phase to a spatial superposition state, which de- 
pends on the momentum of the colliding gas particle and 
is therefore random. This article hence aims at a fun- 
damental change of the understanding of the collisional 
decoherence process. 

Because we are interested in a general understand- 
ing of the collisional decoherence process, rather than 
in details depending on a particular interaction model, 
we use the simplest possible model. That is, we assume 
that the tracer particle as well as the gas particles only 
move in one dimension. Further, the gas is in a thermal 
state using Boltzmann statistics, and we apply the low 
density and high temperature limit. Within this limit, 
each collision is an independent event and we can neglect 
three particle collisions. The gas particles do not interact 
with each other (ideal gas), and their interaction poten- 
tial with the tracer particle is of the hard core type, i.e. 
V(x — x g ) = hnia^oo aS(x — x g ) where the index g labels 
the gas particle. 

Let us briefly review a single collision following refer- 
ence [14 J. The effect of a collision on the tracer particle 
depends on the state of the colliding gas particle. There- 
fore, we start this discussion with a convenient convex 
decomposition of the thermal density operator of a gas 
particle. A particular useful convex decomposition was 
given by Hornberger and Sipe [5 in terms of Gaussian 
minimum uncertainty wave packets \x g ,p g ) o . with 



( 



x n \ x, 



g>i 



>l = 



-iX g Pg/2h 



TT(J n 



e lx' g Pg/h e - 



(Xg 



\fl^l 



(l) 



where x g and p g label the mean position and momentum 
of the wave packet, respectively, and o~ g labels the posi- 
tion uncertainty. It was shown that the density operator 
can be written as 

2nh //v x ,. 
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Here, firiPg 



e- p V 2m ° kBT / ^2irni g k B T is the 



Maxwell-Boltzmann distribution, L is a normalization 



length which is taken to infinity, and T a —T— 2m \ a 2 • 
The reason for a lower temperature in the Maxwell- 
Boltzmann distribution in Eq. (3| is that part of the 
thermal energy of the gas particle has been transferred 
to being a contribution to the momentum uncertainty of 
the states \x g ,p g ) . 

With Eq. ([3| at hand, we can assume that every gas 
particle is in a minimum uncertainty state with position 
uncertainty <j g , while the probability density for a par- 
ticular combination of x g and p g is given by rigfiT^ {Pg), 
where n g is the particle density of the gas. 

For a complete collision of the gas particle wave packet 
with the tracer particle wave packet, it is required that 
the velocity uncertainty of the gas particle state \x gi p g ) a 
is small compared to the relative velocity of the two col- 
liding particles. It was shown in [14] that this is the case 
(at least for most gas particles) if 2m g kBT<j g l ^> h 2 , and 
therefore we will choose o~ g sufficiently large and approx- 
imate HT a (p g ) by /JiriPg) m the following. Furthermore, 
to avoid the discussion of three particle collisions, we re- 
quire o~ g n g <C 1. Therefore, the position uncertainty has 
to satisfy 
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which is generally possible in the high-temperature and 
low- density limit 
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Note that this limit must also be satisfied to use Boltz- 
mann statistics to describe an ideal gas. 

It was further shown in [14J, that, under the additional 
assumption of a slow (compared to the gas particles) 
tracer particle, the collision rate is 



R = n g \/2k B T/^/7mig 



(6) 



It is well understood that if the tracer particle is ini- 
tially (at time t — 0) also in a minimum uncertainty 
state \x,p) a , but with a position uncertainty a related to 
the gas particle's position uncertainty via their relative 
masses according to 
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then a collision results in the remarkable simple product 
state Q31 OS] 



Ug(t)\x g ,p g ) a ®U(t)\x,p) a . 
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Here, t is some time after the collision, and U(t) is the 
free evolution operator of the single particle. The mean 
positions and momenta after the collision relate to the 
initial values according to 
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FIG. 2: (color online) Position (a) and momentum (b) prob- 
ability distributions for the tracer particle during a collision 
in the reference frame of center of mass. While the position 
probability distribution "flows", the momentum distribution 
"jumps" from the initial value to the final one. The solid 
lines in (a) are the corresponding classical trajectories. The 
collision time t c is indicated in (b). 



which are precisely the same relations as in a collision 
of classical particles. The relative mass is denoted by 



m, 



j/m. The final state rt8j) was derived in 



by 



two different methods, and can also easily be confirmed 
by using the scattering operator for a one dimensional 
hard core interaction, i.e. S \p g ) \p) = — \p g ) \p). 

The reader should note that at no time the position 
probability distribution of either particle changes discon- 
tinuously, as might be suggested incorrectly by Eq. ([9| 
and Eq. (10). The actual position probability distribu- 
tion of the state (pi) also includes the free particle evolu- 
tion operator, and in [14] it was shown that the evolution 
of the tracer particle in position representation is com- 
pletely continuous, as can be seen in figure [2] (a). 

We should point out that in this work we neglect the 
possibility of incomplete collisions of the colliding wave 
packets. Therefore, our work is valid on a coarse grained 
time scale 
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where t c is the collision time defined in [14] and shown 
in figure [2] (b). For this very reason, we will use a finite 
time interval (0, t) to study the decoherence process. 



II. REMARKS ABOUT COLLISIONAL 
DECOHERENCE 

Coherences of a quantum state are often described by 
the off-diagonals of the density matrix. This leads to 
the question: What basis should we use to examine de- 
coherence? The two bases which first come to mind are 
the position basis and the momentum basis, and indeed, 
these are the bases usually used in the literature [3]- 
HJ [TTJ [12j [15] . Nevertheless, they are not without prob- 
lems. 

If one uses e.g. the momentum basis, one would be in- 
terested in how long a superposition of the form (\p a ) + 



\pb))/y/2 survives, or, equivalent ly, how fast the coher- 
ences \p a )(Pb\ decrease. There is a major problem in this 
sort of question: A momentum eigenstate (or any state 
which is very localized in momentum) is itself a coher- 
ent superposition of widely spread position eigenstates, 
and one should expect that these position coherences of 
each individual momentum eigenstate do vanish on the 
same time scale, or even faster, as the momentum coher- 
ences of interest! The same reasoning applies to using 
the position basis. 

Our reservations are clearly related to the concept of 
a pointer basis [2j [17]. Pointer basis states should be 
fairly robust to decoherence. That is, if the system is 
prepared in one of these pointer states it will stay there 
for some time, whereas if the system is in a superposition 
of two pointer states, then the coherences typically de- 
crease rapidly in time. These properties makes a pointer 
basis the basis of choice to study decoherence effects. Be- 
cause a momentum eigenstate is highly unrobust to po- 
sition decoherence, the momentum basis is not a good 
pointer basis. Similarly, a position eigenstate is unro- 
bust to momentum decoherence, and therefore also not 
an appropriate pointer state. 

To add some weight to our concern, we have a brief 
look at the decoherence rate found by [11] for an ini- 
tial superposition (\p) + \—p))/y/2 of momentum states, 
i.e. their equation (67). Inserting the definitions (27, 
31, 32) of [11] and assuming that the velocity of the 
tracer particle is small compared to an average gas par- 
ticle, their equation (67) leads to the decoherence rate 

D p = ^ arig ^.fe, where we changed the nota- 
tions according to ours, and a is a constant scattering 
cross section. This formula strongly opposes physical in- 
tuition because it predicts a decrease of the decoherence 
rate upon increasing the temperature, despite the fact 
that an increase of the temperature leads to more pow- 
erful and more frequent collisions. 

Could we possibly single out a pointer basis from a 
measurement interpretation of a single collision? We 



showed in [14] that a colliding gas particle 



v g-> 



P 9 ) an per- 



forms a smeared out measurement on the tracer particle 
in the basis \x,p) a , indicating that these states could be 
used as pointer basis. We know from Eq. p|, that we 
have a choice in the width a g of the gas particle states. 
Therefore, we can use minimum uncertainty states of any 
width a = ^m/m n o- n as pointer basis, as long as a n sat- 
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Q. 

of a single collision. 



v g u g as pointer basis, as long as a g 
isfies relation (|4|), which was required for the treatment 



As is stated in e.g. [17 , there are still many open ques- 
tions about the emergence of a pointer basis from the 
coupling to an environment. Nevertheless, because of 
the heuristic reasoning in the previous paragraph, and 
the lack of sensible alternatives, we will indeed discuss 
decoherence by using superpositions of Gaussian states 
\ x a:Pa) a + \ x b,Pb) a : commonly referred to as cat states 
in the literature. We should note, that [12] superimposed 
several momentum eigenstates to produce a state simi- 



lar to a cat state to study decoherence using quantum 
trajectory theory. This study was limited to the situa- 
tion where the gas particle mass equals the tracer particle 
mass, and each collision leads to a complete loss of deco- 
herence. 

We will need the transformation of an initial (unnor- 
malized) cat state \x a ,p a ) o . + \xb,Pb} a of the tracer par- 
ticle upon a collision with a gas particle in the state 
\ x giPg)a ' ^ using the linearity of quantum mechan- 



ics as well as Eqn. (|8])-(iO), and after tracing out the 
gas particle, we find for the density matrix of the tracer 
particle after a collision 

p(t) = U(t)[\x a ,p a ) a (x a ,p a \+Ce~ % ^\Xa,pa) a (Xb,Pb\ 
+ Ce t(p \x b ,p b ) a (x a ,pa\ + \Xb,Pb) a (x h ,p h \]U\t). 

(12) 



Here, we used x a = x(x g ,x a ) etc. given by Eq. (10), as 
well as 
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which is smaller than one because of a reduction of coher- 
ences due to a measurement performed by the colliding 
gas particle, and 
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is a phase shift induced by the collision. We further de- 
fined the average position xa = (x a -\-Xb)/2 and the posi- 
tion difference xp = (x a — Xb) (analogous for momenta) 
of the two cohering wave packets. 

Note that the phase shift ip = cp(x g ,p g ) depends on 
the state of the colliding gas particle. If this state is 
random as in a thermal gas, the phase shift leads to phase 
averaging and therefore to a loss of coherences. 

Rather than using a density operator representation in 
terms of Gaussian states, we found a Wigner function 
representation (see e.g. [20]) more graphic to display co- 
herences of cat states. Decoherence can then be discussed 
in terms of the vanishing of the oscillatory behavior of 
the Wigner function. In particular, we will compare the 
Wigner function without a collision, to the Wigner func- 
tion with a collision. This way, we will obtain the 'deco- 
herence per collision', which can be multiplied with the 
collision rate Eq. ([6]), to obtain the decoherence rate. 

We wish to mention that a Wigner function description 
of quantum Brownian motion was already put forward 
previously in a very heuristic derivation [18], as well as 
in a more precise approach, but limited to states of the 
tracer particle which are close to a thermal state [19]. 
Both articles are concerned with the general form of a 
partial differential equation for the Wigner function, and 
do not study decoherence. 



We will need the Wigner function corresponding to a 
density operator of the form 

P = \Xa,Pa) a (Xa,Pa\+Ce~ tiP \Xa,Pa) a (xb,Pb\ 

+ <X %(p \Xb,Pb) a (Xa,Pa\ + \Xb,Pb) a (Xb,Pb\' (15) 

Here, c is bound between zero and one and is a measure 
for the strength of the coherences between \x a ,p a ) a and 
\ x b-,Pb)cri an d if determines the relative phase between 
these states. In particular, an incoming tracer particle 
state \x a ,p a ) a + \xb,Pb) a before a collision corresponds 
to c = 1 and cp = 0. 

Using standard methods [20 , the Wigner function cor- 
responding to the density matrix Eq. ( 15 ) is easily found 
to be 
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As is well known, the strength of the coherences, indi- 
cated by oscillatory behavior of the Wigner function, does 
not change due to the unitary free evolution. Therefore, 
in the following, we will mostly use the interaction pic- 
ture [2T] . We can then use the general formula Eq. (16), 
with c = 1 and (p — 0, to plot the Wigner function of 
a state \i/j) = \x a ,p a ) a -\- \xb,Pb) a without a collision in 
figure [3] (a). 

According to Eqs. (T2][l4), a collision with a gas parti- 



cle state \xg,p g ) results in 

( x a/b,Pa/b, x A/D,PA/D) ~^ (x a /b, Pa/b, X A/D, PA/d) 
C = 1 -» C 

p = -+ <p, (17) 

where x a = x(xg,x a ),etc are given by Eq. (TlOJ) , and c 

and (p by Eq. ( 13 ) and Eq. (14), respectively. The Wigner 
function after the collision is plotted in figure [3] (b). 

It is quite astonishing, that, despite choosing a gas par- 
ticle with only four percent of the mass of the tracer parti- 
cle, and, despite using a superposition of very close Gaus- 
sian wave functions, almost all coherences are lost after 
a single collision. If we had separated the initial Gaus- 
sians only slightly more, or had chosen only a slightly 
heavier gas particle, the coherences would be not visi- 
ble at all, because c in Eq. (13) decreases exponentially 
with these parameters. This observation is independent 
of the initial momentum and position of the colliding gas 
particle, as well as whether the Gaussian wave functions 
are separated predominantly in position or momentum. 
It is therefore fair to say, that, unless the gas particle 
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FIG. 3: (color online) The Wigner function of the initial co- 
herent superposition \ip) = \x ai Pa) a + \xb,Pb) a (a), and of the 
state resulting from a collision with \x g ,p g ) a (b). Parame- 
ters are: x a = 15, Xb = 0, p a = 0, pb = 1.5, cr = 4, m = 
I, H=l, x g = 100, £><? = -1, and a = 0.04. 




FIG. 4: (color online) As in figure[3] but with different param- 
eters for the gas particle: p = —0.2, x g = 500, and a = 0.002. 



is much lighter than the tracer particle, the decoherence 
rate equals the collision rate. 

This result becomes even more pronounced, if we av- 
erage over different initial gas particle positions and mo- 
menta, and we will study this effect in the following sec- 
tion for extremely light gas particles, for which the deco- 
herence per collision due to information exchange of the 
colliding particles, (1 — c), will be small. 



resolved. The reason is that the measurement performed 
by the colliding gas particle is so imprecise, that it can 
not distinguish between the two Gaussian wave functions 
of the tracer particle (22] . 

The main effect of a collision with a light gas par- 
ticle is a shift of the entire Wigner function in phase 
space. Of particular interest is that the oscillating part 
just shifts from (xa,Pa) to (xa,Pa), without acquiring 
an additional phase shift (i.e., the relative heights of 
the oscillating peaks do not change). This feature can 
be explained by looking at the argument of the cosine 
in Eq. (16). The sum of the first two terms does not 
change in a collision, because (xaPd— paXd)/(2H) equals 
(zAPD—pA%D)/(2?i')+ ( P, with (p taken from Eq. (14). The 



third and forth term account for a phase shift correspond- 
ing to the shift of the Gaussian, i.e. xa — x f — >> xa — x f 
and p A -p' -^Pa ~p' - 



Of course, if the gas particle state 



iVg, 



is taken from 



a thermal gas, we have to average over all gas particle mo- 
menta p g , weighted by the Maxwell-Boltzmann distribu- 
tion HT{Pg), as well as over all the gas particle positions 
x g which can reach the tracer particle in a given time in- 
terval (0, t). Because each possible combination of x g and 
p g results in a different shift of the Wigner function in 
phase space, it is clear that this procedure strongly sup- 
presses the oscillations. This is the decoherence effect 
we referred to as 'phase averaging' in the introduction. 
It can suppress coherences quickly, even if the measure- 
ments which the gas particles perform are very weak. 

Let us close this preliminary discussion with a note re- 
garding the collision rate. In the limit of a small mass 
ratio a, the tracer particle will be very localized com- 
pared to the gas particles, in both, position and velocity, 
as is evident from Eq. (|7|). Therefore, we do not have 
to use the full rate operator formalism developed in [14] , 
but we can simply assume that the tracer particle is rea- 
sonably localized somewhere near the origin, and a gas 
particle collides with the tracer particle during a time 
interval (0, t) exactly if 



< -XgTTlg/Pg < t 



(18) 



is satisfied. 



III. COLLISIONAL DECOHERENCE FOR 
LIGHT GAS PARTICLES 

In this section, we discuss the more interesting situa- 
tion, when a collision only partially destroys the coher- 
ences, which is the case if the colliding gas particle is 
extremely light (mass ratios much smaller than one per- 
cent, as we will see). Then, we need a more quantitative 
measure of the decoherence process, which we will de- 
velop below. 

Figure [4] shows how the Wigner function of a super- 
position state changes due to a collision with a light gas 
particle. The coherences after the collision are still well 



A. Position decoherence 

To study position decoherence, we consider an initial 
tracer particle state \x a ,p} a + \%b,p} (T ' An example of the 
corresponding Wigner function is plotted in figure [5] (a). 
The Wigner function after a collision with a gas particle 
state \xg,Pg) a is obtained following section III) but to 
study the decoherence due to a thermal gas, we average 
over different initial positions x g and momenta p g of the 
gas particle. 

The Wigner function W Pl after a collision, averaged 
over 200 pairs (x g ,p g ) taken from an appropriate distri- 
bution (see below), is shown in figure [51(c), (e), and (g), 



Decoherence Rate 




we can use c 



1, and concentrate on the cosine of the 



FIG. 5: (color online) (a): The Wigner function for the tracer 
particle before a collision, (c), (e), (g): The average Wigner 
function after one collision with a gas particle at temperature 
T = 0.2, 0.5, and 1.5, respectively, (d), (f), (h): The change 
of the Wigner function due to a collision at temperature T = 
0.2, 0.5, and 1.5, respectively, (b): The relative change of the 
Wigner function at the origin due to a collision. This serves as 
a quantitative measure of the 'decoherence per collision'. The 
solid line is the first order expansion in temperature, Eq. (27). 
Parameters are: x a = 20, Xb = —20, p a = 0, pb = 0, a = 
4, m = 1, h = 1, k B = 1, t = 20, and a = 0.0001. 



for the respective temperatures T = 0.2, 0.5, and 1.5. 
The corresponding changes of the Wigner function from 
the initial one (figure [5] (a)) are shown figure [5] (d) , (f), 
and (g). 

As a quantitative measure of coherence, we use the 
hight of the maximum peak of the oscillations. In partic- 
ular, the relative change of the maximum peak serves as 
'decoherence per collision', and is plotted over the tem- 
perature T in figure [5] (b) . 

In the following, we derive an expression for the 'deco- 
herence per collision'. In the limit of a light gas particle 



Wigner function Eq. ( 16 ) 
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where we used pjj = and Pa = P for our choice of the 
initial cat state. 

As noted earlier, a collision does not change the sum 
ip — pxd/(2K). Further, we use 1 + a « 1 in the limit 
of light gas particles, which leads to xjj ~ xjj and p ~ 
p + 2p g . Therefore, we find the oscillating term of the 
Wigner function W Pl after a collision by averaging over 



cos 
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2h 
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(20) 



In particular, we use the reduction of the maximum of 
the oscillations as a quantitative measure of decoherence. 
Because this maximum before a collision is at p' given by 
cp—px]j/(2h) +xjj(p — p')/h = 0, we substitute this into 
Eq. ( [20| ), and the coherences after one collision are then 
obtained by averaging over 



cos 



2X D Pg 



(21) 



We see that the relative phase 2p g xo/h added by the 
collision depends solely on the momentum of the colliding 
gas particle, and we will need the momentum probability 
distribution of the colliding gas particles to perform the 
averaging over Eq. ( [21] ). 

According to Eq. pF, the probability density of finding 
a gas particle in the state \x g ,p g ) is given by n g HT{Pg), 
where n g is the particle density of the gas. Knowing that 
a gas particle of momentum p g collides with the tracer 
particle exactly if the position x g satisfies Eq. ( fl8] ), we 
can write down the normalized momentum probability 
distribution of a colliding gas particle 



C( Pg ) = 



\Pg\ 



2m q kBT 



exp 



2m g kBT 



(22) 



Note that this distribution does not depend on the length 
of the considered time interval (0, t). This will lead to a 
'decoherence per collision' (and therefore to a decoher- 
ence rate) which is independent of the considered time 
interval, as should be expected in a Markovian process. 
We finally find for the coherences after one collision 
with a thermal gas particle 



cos 



2x D p g 



C( Pg ) 
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Jo 



d Pg 
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• aknT 



exp 



2m q k B T 



cos 
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(23) 
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ft 
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sin ' 



2m q k B T 



ft 



(24) 



where we used integration by parts. We deduce for the 'decoherence per collision' of position superposition states 

(25) 



Collision 
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FIG. 6: (color online) (a): The 'decoherence per collision' 
plotted over the temperature, for parameters as in figure [5] 
but for larger temperatures. The dots are from the numerical 
Wigner function at the origin, and the line is from Eq. ( |24| ). 
The small discrepancy is due to the neglected spreading of 
the Gaussians due to the averaging process. 
(b): The Wigner function after a collision at temperature 
T — 7 shows indeed oscillations in opposite phase to the initial 
Wigner function (see figure [5] (a)). This leads to a decoherence 
rate which exceeds the collision rate. 



An interesting limit to consider is 2x B ^/2rn g ~kBT <C ft, 
in which the 'decoherence per collision' is small. This sit- 
uation corresponds either to a small position separation 
xd, or to a small average momentum transfer y^m^/c^T, 
and was studied in the seminal work of Joos and Zeh 3 . 



The integrals in Eq. (25)- (26) can then be carried out 
and we find 
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Collision 



4m Q kBT 9 

— — x 
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Sng^/ffTg-jkBT) 3 / 2 

^h 2 l 



(27) 
(28) 



The decoherence rate D x agrees up to some constant with 
the one found by [3 for three dimensional collisions. 



The 'decoherence per collision' Eq. (25) is plotted over 
temperature in figure [6] (a) for the same parameters as in 
figure [51(b), but for higher temperatures. It might come 
as a surprise, that the decoherence rate exceeds the colli- 
sion rate for xr>\/2m g kBT '/ft > 1. The reason is that in 
this regime, the Wigner function after a collision shows 
oscillations, which are out of phase with the oscillations 



of the initial Wigner function, as shown in figure [6] (b). 
Therefore, if we write down the actual (interaction pic- 
ture) Wigner function after some small time t as 



W, 



Pit) 



(1 - Rt)W P0 + RtW p 



(29) 



where R is the collision rate Eq. ([6]), and po and p\ are the 
density operators corresponding to no collision and to one 
collision, respectively, then the oscillations in W Po and 
W Pl interfere destructively. These out-of-phase coher- 
ences in turn can be understood by noting that the mo- 
mentum distribution of the colliding particles, Eq. (22), 
is not peaked at p g = 0, but rather at p g = ±y / 2m^T. 

Having found the 'decoherence per collision' due to 
phase averaging, we can draw a quantitative comparison 
with the 'decoherence per collision' due to information 
exchange, which is 1 — c « ax 2 D /cr 2 . Because of the first 
inequality of Q, this decoherence effect is indeed negli- 
gible if the density and temperature of the gas are such, 
that it can be considered an ideal Boltzmann gas. This 
is also true in the regime discussed in section [TTJ because 
phase averaging is sufficient to remove any coherences in 
a single collision. 

We note that it is often stated in the literature 
[3j [15] that, if the separation xd of two interfering 
wave packets is larger than the thermal wave length 
A = ft/ y/27rm g kBT of the gas, then a colliding gas par- 
ticle can distinguish between the two interfering wave 
packets, therefore removing their coherences. In finding 
that the decoherence rate is about the collision rate if 
x d ^ V^V we confirm the latter part of this statement, 
but we also show that the loss of coherence is by no means 
related to a measurement performed by the gas particle, 
but due to classical phase averaging resulting from the 
randomness of the momentum transfer. 



B. Momentum decoherence 

In this subsection, we will show that the decoherence 
of momentum superposition states is not a direct process. 
Instead, two coherent wave packets with momentum sep- 
aration pd will, after some time, acquire a position sep- 




FIG. 7: (color online) (a): The Wigner function for the tracer 
particle before a collision, (c), (e), (g): The average Wigner 
function after one collision with a gas particle within a time 
interval (0,£) at temperature T, with t = 20, T — 0.5 in (c); 
t = 50, T = 0.5 in (e); t = 20, T = 3 in (g). (d), (f), (h): 
The corresponding change of the Wigner function due to a 
collision, (b): The relative change of the Wigner function at 
the origin due to a collision as function of the time interval 
for T = 0.5 as well as T = 1. Parameters are: x a =0, Xb = 
0, p a = 1.2, pb — —1.2, a — 4, m — 1, h = 1, /c^ = 1, and 
a = 0.0001. 



aration #£> = p^t/m, which then leads to position de- 
coherence. Any direct momentum decoherence will turn 
out to be negligible in a high temperature and low density 
gas. 

We consider the initial tracer particle state \x,pa) . + 
\ x iPb} a , whose Wigner function is plotted in figure ul (a). 
Again, the main source of decoherence will be phase av- 
eraging. Contrary to the previous subsection, where the 
relative phase of the two Gaussian wave packets after a 
collision depended on the initial momentum of the collid- 



ing gas particle, for momentum decoherence this phase 
depends on the initial position of the colliding gas parti- 
cle [23] . As the variation of the initial position of a collid- 
ing gas particle is increased, the more 'decoherence per 
collision' we will find. For a given gas particle momenta 
p g , the gas particle position can be anywhere within the 
interval (—p g t/m g ,0). Therefore, the 'decoherence per 
collision' will not only depend on the temperature T (for 
the distribution of p g ) : but also on the considered time 
interval. In figure [T] (c) - (h), the effects of a collision on 
the Wigner function is shown for different temperatures 
and time intervals. The dependence of the 'decoherence 
per collision' for low temperatures and short times turns 
out to be linear in temperature and quadratic in time, as 
shown in figure [71(b). As a result, it is not possible to de- 
fine a time independent decoherence rate for momentum 
superpositions. We will provide a physical interpretation 
at the end of this subsection, and first give a mathemat- 
ical explanation of these results. 

For this purpose, we consider again the oscillating co- 
sine within the Wigner function Eq. (16) 



<P- 



2h 



x — X 



PD- 



h 



(30) 



at its maximum cp + xpd/(2H) — Pd(x — x')/h = 0. As 
discussed before, a collision does not change the sum 
if + xpd/(2K) in the cosine, and we only have to con- 
sider the change of x to x + 2ax g in the last term of the 
cosine. Therefore, we find the 'decoherence per collision' 
by averaging over cos(2ax g pjj/h). Here, we need the 
normalized probability distribution C(x g ) of the initial 
position of the colliding gas particle, which is obtained 
from the probability density n g HT(Pg) by integration over 
all p g , for which a gas particle with position x g can reach 
the tracer particle (i.e. 
upper sign is for positive x g ) 



x g ^ Pgt/rrig ^ =F oo , where the 



C(x g ) cx / dp g exp\ r-— - . 

y/27rm g k B T J\x g \m g /t \ 2rn g k B T J 

The distribution is normalized either by integration over 
x g , or directly by dividing by the collision probability Rt. 
After substituting u = p g j \j2m g kBT we find 

due~ u \ (31) 



C(Xg) 

and therefore 



ty/2k^T J\^A 



ty/2k B 3 



(2ap D x g 
cos — - 

V ft 



C(x g ) 



t^JW^F 



dx n 



2ap D x g 

h 



\ x g\\ 



due 



t^2k B T 



m h 



tyj2m g k B Tp D Jo 



» _ n 2 . ( 2ty/2m g k B Tp D \ 

due sin — u . 

\ ma / 



(32) 



Because this function represents the coherences after one collision, we have to subtract it from unity to obtain the 
'decoherence per collision' 
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(33) 
(34) 
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The solid lines in figure [7] (b) are taken from Eq. (33), 
and agree well with the data (dots) obtained from the 
numerical Wigner functions directly. The approximation 
in Eq. (34) is valid if the 'decoherence per collision' is 
small. 



Similar to position decoherence, we see from Eq. (34) 
that momentum decoherence due to information ex- 
change (1 — c « a(j 2 pjj/h 2 ) is negligible (unless for tem- 
peratures and times so small, that relation ( pTj ) is vio- 
lated). Hence, we established that also momentum deco- 
herence is due to phase averaging. 

At first, the increase of the decoherence rate with 
the considered time interval seems to be at odds 
with the uniformity in time in the following sense: 
If we split a time interval (0, t) into sub intervals 
(0, t/N), (t/N, 2t/N), ...,(£- t/N, i), the decoherence 
rate of the entire interval should be the averaged deco- 
herence rate of all the sub intervals. If we now assume 
that the decoherence rate for each sub interval is the same 
("uniformity in time"), we would be lead to the conclu- 
sion that the decoherence rate for the interval (0, i) equals 
the decoherence rate for the subinterval i0, t/N), clearly 
contradicting Eq. (33) and Eq. (34). 



In the above argument, we made the following concep- 
tual error: by assuming the same decoherence rate for 
each sub interval, we implied the same initial cat state at 
the beginning of each sub interval. But instead, by the 
time the n-th sub interval starts, the cat state evolved to 
U(nt/N)(\x,p a ) a + \x,pi ) ) a ) : and acquired a position sep- 
aration x B = point /N)/m. Therefore, we have to add 
position decoherence for all but the first sub intervals. 

In the limit N — >• oo, there is no momentum deco- 
herence in the infinitely small sub intervals at all, but 
instead, a continuously increasing position decoherence 
according to Eq. (25), with #£>(£') = p B t f /m. Indeed, 
substituting this time dependent position separation into 
Eq. (25) and averaging over all times t' E (0,t), one ex- 
actly recovers Eq. (33) [24] . 

In other words, the decoherence which an initial cat 
state \x,p a ) + \x,pb) experiences during a time interval 
(0, t) is perfectly explained by position decoherence of the 
evolving state U(t f )(\x,p a )-\-\x,pb)). This leads us to the 
physical interpretation that momentum decoherence is 
not a direct process, but results indirectly from position 



decoherence due to position separation which the tracer 
particle acquires over time. 



IV. CONCLUSIONS FROM THE STUDY OF 
DECOHERENCE 

We showed previously [14] that in one dimensional col- 
lisional decoherence, the quantity measured by a colliding 
gas particle does not only depend on physical parameters 
like density and temperature, but also on the choice of de- 
composition Eq. ([3| of the density operator of a thermal 
gas particle. It is therefore reassuring to find in this arti- 
cle, that measurement effects as a source of decoherence 
are negligible in the high temperature and low density 
limit, where Eq. ^ is valid. The reason is that mea- 
surement effects are small compared to phase averaging 
effects, which arise from a random relative phase added 
to a superposition state during the collision process. 

We further arrive at a neat interpretation of the deco- 
herence process of a superposition of two Gaussians wave 
packets. The decoherence due to a collision depends on 
the position separation of the two Gaussian wave pack- 
ets at the time of the collision. In contrast, there is no 
direct decoherence due to the momentum separation p B 
of the two coherent wave packets. Instead, over time, 
the momentum separation changes the position separa- 
tion according to #£>(£) = x B + tp B /m. This leads to an 
indirect influence of p B on the decoherence rate, which, 
if there is no initial position separation, is described by 
Eq. (pi). 



Further work is required to see whether this drastic 
change in the understanding of the collisional decoher- 
ence process also applies to three dimensional systems. 
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